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Abstract 

In this paper we investigate gravitational instability of shocked gas layers using linear analysis. An 
unperturbed state is a self-gravitating isothermal layer which grows with time by the accretion of gas 
through shock fronts due to a cloud-cloud collision. Since the unperturbed state is not static, and cannot 
be described by a self-similar solution, we numerically solved the perturbation equations and directly 
integrated them over time. We took account of the distribution of physical quantities across the thickness. 
Linearized Rankinc-Hugoniot relations were imposed at shock fronts as boundary conditions. The following 
results are found from our unsteady linear analysis: the perturbation initially evolves in oscillatory mode, 
and begins to grow at a certain epoch. The wavenumber of the fastest growing mode is given by fc = 
2^/2'kGpe.M/cs, where pE, Cs and M are the density of parent clouds, the sound velocity and the Mach 
number of the collision velocity, respectively. For this mode, the transition epoch from oscillatory to 
growing mode is given by tg = 1.2/^/2ttGpbM- The epoch at which the fastest growing mode becomes non- 
linear is given by 2ASo"-^/\/2ttGpeM, where Sq is the initial amplitude of the perturbation of the column 
density. As an application of our linear analysis, we investigate criteria for collision-induced fragmentation. 
Collision-induced fragmentation will occur only when parent clouds are cold, or uq = ^c^R/2GM < 1, where 
R and AI are the radius and the mass of parent clouds, respectively. 
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1. Introduction 

Shock waves that propagate in interstellar gases play 
important roles in the formation of stars and other struc- 
tures in the universe. Shock waves arc formed by cloud- 
cloud collisions (Hasegawa et al. 1994, Sato et al. 2000), 
supernova explosions (Herbst & Assousa 1977), expan- 
sions of HII regions by ionizing photons emitted from 
OB stars (Blaauw 1964), and so on. Shock waves sweep 
out external gases and form highly compressed sheet-like 
structures behind the shock fronts. These layers, which 
are sheet-like structures, are often supposed to fragment 
and eventually to form stars. 

Shock compression contributes especially on the for- 
mation of aggregations of stars. Elmegreen & Lada 
(1977) proposed a scenario of a sequential formation of 
OB star subgroups. The expansion of a HII region that 
is formed by an OB subgroup sweeps up the gases. A 
highly compressed layer forms, and then possibly frag- 
ments. Thereafter, a new OB star subgroup forms from 
the layer, and the ionizing photons emitted from the new 
subgroup induce the formation of next subgroup by the 
same process. Kumai et al. (1993) proposed a scenario of 
a globular cluster formation that is induced by the gravita- 
tional instability of a layer formed by a supersonic collision 
between large proto-grobular cluster clouds. Observations 
showed that galaxies containing young star clusters have 
larger random velocity dispersions of gases than those that 
do not have young clusters (Fujimoto & Noguchi 1990). If 
the large random velocity dispersion of gases is stimulated 



by mergers or interactions between galaxies (Ashman & 
Zepf 1992), the scenario of star cluster formation by 
cloud-cloud collision is indicated. 

In the process of star cluster formation from the layer, it 
is crucial to understand when and how a shock-compressed 
layer fragments. Goldreich & Lynden-BcU (1965) (here- 
after GL65) and Elmegreen & Elmegreen (1978) investi- 
gated the stability of a static isothermal self-gravitating 
layer that is supported by constant external pressure on 
both sides using linear analysis. They showed that the 
wave length of the most unstable mode is approximately 
the scale height of gravity, Hq = Cg/ y/2TTGpQo when the 
density at the central plane, poOi is larger than the bound- 
aries, pb; of the layer (poo/Ph > 1)- When poo/Pb — 1, the 
wave length of the most unstable mode is the thickness of 
the layer. The maximum growth rate is about the free-fall 
growth rate of the layer in both cases. 

However, a layer formed by shock compression is 
bounded by shock waves on at least one side. The effects 
of the shock boundary strongly influence the stability of 
the layer when poo/Pb — 1- This is because the layer is 
mainly confined not by self-gravity, but by external pres- 
sure. Vishniac (1994) (hereafter V94) derived a simple 
dispersion relation of a layer bounded by shocks on both 
sides using shock boundary conditions, while neglecting 
the time dependence of the column density. The disper- 
sion relation of the shocked layer is quite different from 
that of the static layer when poo/pb — 1 (see Appendix 
2 and 3). Therefore, the results with a constant pressure 
boundary are not applicable to a layer bounded by a shock 
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Fig. 1. Schematic picture of a compressed layer between the 
shock fronts. 



front when poo/Pb — 1- Since the layer experiences the 
state of Poo /Pb — 1 in early phase from the moment when 
the layer is formed by collision, an analysis that takes ac- 
count of the effect of the shock boundary is important. 

In the analysis of V94, the time evolution of the column 
density is not considered. In this paper, in order to con- 
sider the fastest growing mode in the layer that is formed 
by cloud-cloud collision, we take into account both shock 
boundary and time-dependent column density. The dis- 
tribution of the physical variables across the thickness is 
also taken into account accurately. 

In section 2, wc set up a model of a layer. Before our lin- 
ear analysis, previous works (V94; Whitworth et al. 1994) 
using a 1-zone model is reviewed. In section 3, an un- 
perturbed state is described. In section 4, perturbation 
equations are derived and a numerical method, boundary 
conditions and initial conditions are shown. The results 
arc shown in section 5. Criteria for collision-induced frag- 
mentation arc discussed in section 6. Our study is sum- 
marized in section 7. 

2. Model 

2.1. Assumptions 

We consider a head-on collision of two physically iden- 
tical clouds with sufficiently large Mach number {A4 ^1). 
The equation of state is assumed to be isothermal. In this 
paper, Cartesian coordinates {x,y,z) are used and the z- 
axis is perpendicular to a layer. The layer is assumed to 
extend infinitely in the x — y plane. We can discuss the sys- 
tem in the x — z plane without any loss of generality in the 
linear regime. Physically identical clouds are assumed to 
collide along the z-axis with velocity ±Ve and density pE 
at t = and z = 0. The two shock fronts at z = ±Zsh prop- 
agate in both directions and a highly compressed layer 
is formed between two shock fronts (see figure 1). The 
isothermal shock boundary conditions are given by 



Psh = PeM , ish 
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and Ve 



M 



1 

M 



(1) 



where A4 = T4h/cs is the Mach number and Vsh is the ve- 
locity of the clouds in the rest frame of the shock front (see 
Appendix 1 in detail). The unperturbed column density, 
Eq, is defined by 



Po{z)dz = 2pEMc^t, 



(2) 



which increases as ex < by gas accretion. Hereafter, we 
investigate gravitational instability of the layer between 
two shock fronts. 

2.2. Dispersion Relation under 1-zone Model 

As an essential guide for our linear analysis, a simple 
dispersion relation using a 1-zonc model is summarized 
(see V94 in detail). In the 1-zone model, we consider av- 
eraged physical quantities across the thickness. The layer 
is assumed to be geometrically thin, i.e., Zsh/Ho <C 1 and 
kLo ^ 1, where Lq and Hq = Cg/ ^finGpsh are the thickness 
and the scale height of the layer, respectively. Averaging 
physical variables across the thickness and neglecting the 
time evolution of Eo(t) with shock boundary condition, 
V94 derived an approximate dispersion relation which is 
given by 

w2-c2fc2_27rGfcI]o(t). (3) 

The derivation of equation (3) is shown in Appendix 2. 
The dispersion relation is identical to that of an infinitely 
thin disk without boundary condition. In Appendix 3, 
we derive a dispersion relation of a static layer which is 
confined by constant external pressure by the using 1-zone 
model for comparison. 

According to equation (3), the time evolution of the 
perturbations with wavenumbcr k is described as follows: 
in the early phase i ^ 0, the equation (3) is approximately 
^ Cgfc^, because So(i) is small. Thus, the perturba- 
tions arc initially in an oscillatory mode. Since c^fc^ is 
constant with time, while 2'!iGkY,Q{t) increases with time, 
the perturbations will change from the oscillatory mode 
to a growing mode at a certain epoch. 

This transition epoch can be derived by a simple or- 
der estimation. The gravitational instability sets in if 
> l/t (Whitworth et al. 1994), where l/t is the 
evolution rate of the unperturbed layer by gas accretion. 
The transition epoch, tc,k, at which perturbations with 
wavenumbcr k begin to grow, is given by 



tc,k 



2^GfcI]o(tc,fe)-c2fc2. 



(4) 



Using the characteristic timescale and the wavenumbcr 
defined by 



tc = 1/ \/2t:GpeM and fee = \/2t:GpeM/cs 



(5) 



equation (4) is rewritten as an independent form on Mach 
number as 



■J— = \/2fc.^^K^ 



(6) 



where t — tjtc and k = kjk^. Figure 2a shows the depen- 
dence of fc,K on the wavenumbcr k. Figure 2a indicates 
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Fig. 2. (a)Dependence of the transition epoch, Tc^k, from oscillatory to growing mode on wavenumbers k in the 1-zone model. The 
curve shows a plot of equation (6). (b)Dispersion relation [equation (3)] at f=l, 2 and 3, where lu = t^u). The solid lines satisfy 
\J —u)'^ > l/t> where perturbation grows faster than growth of the column density. The open circle indicates the mode that begins to 
grow at the earliest epoch, k = 1. 



that the mode with k = 1 begins to grow at the earliest 
epoch f = 1. Then, the perturbations with k ~ 1 switch 
to the gravitational instability successively. Whitworth 
et al. (1994) estimated fragmentation mass scale using 
K = 1. However, the earliest transition mode, k = 1, is 
not the fastest growing mode in later epoch. Figure 2b 
represents the dispersion relation [equation (3)] at f =1, 2 
and 3, where Co = t^iLi. The open circle indicates the ear- 
liest transition mode, «= 1. The gravitationally unstable 
state that satisfies \J —uP- > 1/i is shown by the solid lines 
at each epoch. Figure 2b indicates that u) and n of the 
fastest growing mode varies with time. 

In the discussions so far, have been based on equation 
(3), which was derived by neglecting the time dependence 
of So(i). However, especially during i<ic,fc, the evolution 
rate of the layer 1/t is larger than \/ —uj^ . Therefore, the 
evolution of So(t) can not be neglected. In this paper, in 
order to accurately predict the fragmentation scale of the 
evolving layer at the earliest epoch, the evolution of the 
unperturbed layer is taken into account consistently. In 
the following sections, we investigate the gravitational in- 
stability in the linear regime with a time-dependent So(i)i 
as well as shock boundary conditions without the 1-zonc 
approximation. 

3. Unperturbed State 

In our linear analysis, the unperturbed state is the time- 
evolving layer formed by cloud collision. To begin with, we 
define the distribution of the density, p{z,t). The density 
distribution of a hydrostatic isothermal self-gravitating 
layer that is supported by external pressure, Pb, is an- 
alytically written as 



Poo 



-, Ho 



cosh^(z/ifo) ' " V27rGpoo 
The relationship among poo; Pb and Sg is given by 



(7) 



Poo 



7rGS| 
2c2 ■ 



(8) 



Equation (8) shows that a set of external pressure Pb and 
the column density Sq provides a unique density distribu- 
tion p{z). 

The sound-crossing time across the thickness Lq and 
evolution time are defined by tcross = Lq/cs and tcvo = 
S]o/(dS/dt) =t, respectively. Since the thickness Lq varies 
as ^ 2cst/A4, the ratio between the two timescales is 
Across /ievo ^ 2/M <?C 1. Therefore, we assume that the 
layer evolves while maintaining the equilibrium density 
distribution in the z-direction with a large Mach number. 

Next, we assume that the Mach number, M at the 
shock front is constant. We consider the evolution of the 
shock front in the rest frame of the z = plane. According 
to equation (1), if Ve is constant, the shock front velocity, 
Zsh, is constant without gravity. If gravity is taken into 
account, as the effect of gravity becomes important, igh 
decreases with time. Consequently, in the rest frame of 
shock front, the velocity of accreting gas, V^h = Ve + ish, 
and the Mach number A4 decrease. However the frac- 
tional change of Mach number is as small as <C 1. 
Therefore, the constant M assumption well approximates 
a constant Ve- 

The velocity distribution in the postshock layer is as- 
sumed to be zero. In reality, Vz{z,t) is not zero in order 
to attain the hydrostatic state at each instant of time. 
However, this assumption is also valid when the Mach 
number is sufficiently large. The validity of this assump- 
tion is discussed in Appendix 4. 

The above assumptions completely provide the time 
evolution of the density distribution and the position of 
the shock front analytically. With a constant Mach num- 
ber, the column density and the external pressure are 
given by 



Eo(t) = 2pEXc,i, Pb^c^pEX^ 



(9) 



Substituting equation (9) into equation (8), the density at 
the center is obtained as 



Poo = /OEX^(l + r2), 



where r is defined by r = ty/2T:GpE, which is the time 
normahzed by the free-faU timescale in the preshock re- 
gion. Equation (10) shows that the layer has approxi- 
mately uniform density distribution, poo/Pb ^1, at r ^ 1, 
and poo/Pb = 2 at r = 1. Because the density at the shock 
front is given by p^M^, the time evolution of Zsh{t) is 
written with equation (7) as 
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sh = Hq cosh ^ VT+T^ 



1 cosh ^\/1 + t2 



M ^2ttGpe VTT^ ' ^^^^ 
When r ^ 1, equation (11) becomes Zsh ~ Cgt/Ai. As t 
increases, the increasing rate of Zsh is suppressed by self- 
gravity. At T ~ 1.5, Zsh reaches the maximum and begins 
to decrease. 

4. Perturbed State 

In this section, we consider perturbations on the unper- 
turbed state, which is defined in the previous section. To 
begin, the basic equations of ideal self-gravitating isother- 
mal fluid are the mass conservation. 



| + v.(H = o, 

the momentum conservation, 
d 

— (pu) + V ■{pv(g)v + PI) + pV(j) = 0, 

the Poisson equation, 

V^-4> = 47rGp, 
and isothermal condition, 

P 



: const. 



(13) 



(14) 



(15) 



where / is the unit matrix and <f> is the gravitational po- 
tential. The perturbation variables {5p, Svz, 6vx and S(j)) 
are defined by 

p{x,z,t) = poiz,t) + Sp{z,t)e''''', (16) 

v,{x,z,t)=v,^Q{z,t) + Sv,{z,t)e''''', (17) 

v,{x,z,t) = 6v,{z,t)e'''^ , (18) 



and 



(bix,z,t) = Mz,t) + 5ct>{z,t)e'''\ (19) 

where we take the Fourier mode in the x-direction. 
Subscript "0" denotes quantities of the unperturbed state. 
By linearizing the basic equations about the perturbation 
variables, the perturbation equations are given by 

dSp d 
dt dz 

9 , - , d 



.^P)^-^P—-Po- 



dz 



dz 



(20) 
(21) 



{po5vx) = ~ikcl5p - ikpoScj), 



- AirGSp, 



(22) 



(23) 



dz^ 

where we assume Vz,o = 0. 

The rippling of the shock front 6zsh is defined by 

Zsh{x,t)^Zsh,Q{t) + Szsh{t)e'''''. (24) 

The column density perturbation, which is defined by 



ST, — I p(x,z,t)dz — Sqi 
can be divided into two contributions, as follows: 



(25) 



= pshSL + 5a, (26) 

where SL = 26zsh and Sa = J^''^'° ^Spdz. The first term 
describes the fluctuation of the shock boundary, and the 
second term describes the contribution of density pertur- 
bation inside the layer. 

Since the unperturbed state is neither static nor self- 
similar solution, the perturbation equations are not re- 
duced to some eigenvalue problems. Therefore, we directly 
integrate the perturbation equations over both time and 
the z-direction to investigate the time evolution of the 
perturbations. 



(12) 4.1. Numerical Method 



The upwind finite-difference method is used as numer- 
ical scheme. Because the perturbation equations have 
source terms, the numerical flux is modified appropriately 
(Hubbard & Garcia-Navarro 2000). 

We performed three tests of our code. First, we calcu- 
lated fc = perturbation without self-gravity. The results 
of our calculation were compared to the results of a one- 
dimcnsional, non-linear hydrodynamical calculation using 
2nd-order Godunov code in the same situation. Our lin- 
ear results agree with the non-linear calculation very well. 
Secondly, we calculated the case with fc ^ without self- 
gravity. V94 derived the dispersion relation of the even 
mode to be w = ±Csfc(l — We obtained the same 

frequency as V94. Thirdly, we calculated stability of the 
static self-gravitating layer. Our results agree with the 
growth rate of GL65. In these test calculations, our code 
provided correct results. 

4-. 2. Boundary Conditions 

The boundary conditions are set at z = and z = Zsh,o- 
We investigate the even mode, and the boundary condi- 
tions at z = are given by 

(t)„o^°' (^)^^„^0,>»-M.^0)^0,.7, 

At the shock front 2sh,0j the physical variables must sat- 
isfy the Rankine-Hugoniot relations. Thus, we impose lin- 
earized Rankine-Hugoniot relations at Zsh,o (see Appendix 
1 in detail), which are given by 



Sps 



Sv, 



dpo 
dz 
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,sh ■ 



dvz 
dz 



dt 
1 



(28) 



dSzsh 



1 + T7I (29) 



dt 



and 



-ikVESzg 



'sh 



dz 

where subscript 
front. 



kS(t>sh + 47rG(psh - PE)5zsh = 0, 



(30) 



(31) 



"sh" indicates the value at the shock 
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Fig. 3. Time evolution of (5S for ft = 4 and M = 10. 
line corresponds to initial phases ^ = n/A, ir/2, 2tt/3, 
respectively. 
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4-. 3. Initial Conditions 

The initial epoch is set to be Tini = 5 x 10"'^, but the 
results are not sensitive to ri„i as long as it is sufficiently 
small. As an initial condition, an oscillatory mode (see 
subsection 2.2) without self-gravity is adopted. The time 
evolution of perturbations depends on an initial phase ip 
of perturbation in oscillatory mode. In figure 3, the evolu- 
tion of SH with different initial phases is shown for A^=10 
and K = 4. In this paper, we focus on the evolution with 
the initial phase by which grows fastest. 

5. Results 

5.1. The Evolution of Perturbations 

Figure 4 shows the time evolution of i5E for various 
Mach numbers, (M = 10, 30, 50 and 100). Figures 4a 
and 4b correspond to k = 1 and 2, respectively. The thin 
dot-dashed line shows the result without self-gravity for 
A4 = 10. The thick solid line shows the result with self- 
gravity for A^=10. In both of figures 4a and 4b, the lines 
for M. = 10 with and without self-gravity are identical ini- 
tially (r < 0.1). Later, the thick solid line with self-gravity 
changes to a growing mode at r ~ 0.2 — 0.3. We define the 



epoch of transition from an oscillatory mode to a grow- 
ing mode as Tg, which corresponds to the inflection point 
of SH (d2,5E/dr2 = 0). We will find that Tg is well ap- 
proximated by ic,fc in the 1-zone model in subsection 2.2. 
In figures 4a and 4b, it is also seen that the growth rate 
is larger and Tg is smaller for the larger Mach number. 
The dependences of the growth rate and Tg on the Mach 
number are discussed quantitatively in subsection 5.2. 

As shown in equation (26), ST, consists of Psh^L and 
5<T. In figure 5, it can be clearly seen that the fluctua- 
tion of the boundary Psh^L dominates in the evolution of 
the perturbations during t ^ 1. On the other hand, S<t 
stays much smaller than PahSL in t ^ 1. This means that 
the layer behaves like an incompressible fluid for t ^ 1. 
For this reason, as an indicator of the perturbation am- 
plitude, we mainly use the column density perturbation, 
SYi, instead of Sp(z). 

Figure 6 shows the density perturbation, Sp{z), in- 
side the layer at various epochs, t = 0.2, 0.4, 0.6 and 
1.0 for Ai = 10 and k = 1. As normalization, we set 
Szsh/zsh,o = 0.05 at T = 0.6. The abscissa and the or- 
dinate axes are normalized by zq = Cs/{\/2t:GpeM) and 
Psh = -M^pe, respectively. The position of the shock front, 
Zsh.o, is seen as the right edge of each line. At t ^ 0.3 the 
perturbations become gravitationally unstable (c.f., figure 
4a), and 5p increases with time. 

Figure 7 shows cross sections of the perturbed layer at 
(a)T = 0.6 and (b)T = 1.0. The normalization of the per- 
turbation is the same as in figure 6. It can be seen that 
the perturbation grows gravitationally from t = 0.6 to 1.0. 
Near the equatorial plane, inflow in the a;-direction is gen- 
erated by gravity, which reflects to enhancing the shock 
front. In contrast, near the shock boundary, the perturbed 
velocity field is outward in the x-direction. This tangen- 
tial flow is the result of the shock boundary condition. 
This momentum flow acts against the force of gravity (see 
Appendix 2.2 in detail). 

5.2. The Growth Rate and Transition Epoch 

We discuss the dependence of the growth rate and the 
transition epoch, Tg, on Mach number. In the 1-zone 
model, according to equation (5), the dependence of the 
transition epoch, Tc^k, on Mach number is A4~^^'^. To 
see whether the same scaling relation holds or not in the 
model without the 1-zone approximation, we replot figure 
4 in figure 8 using f = rAi^^^ as the abscissa. 

In figures 8a and 8b, it is seen that each line for different 
Mach numbers is nearly identical. This result is explained 
as follows: in the 1-zone model, from equation (3), the 
time evolution of 5Y, is expressed as 



dlnSH 
dt 

where f = 
obtain 



dlnSH 
df 



V2h 



(32) 



A4t. Integrating this equation over f, we 



SH cx exp 



3k 



|3/2 



(33) 



Equation (33) shows that the time evolution of SH seems 




Fig. 4. Time evolution of <5E for the case with (a)K = 1 and (b)K = 2. In each panel, thick lines correspond to the case with 
self-gravity {Ai=10, 30, 50 and 100), while the thin dot-dashed line corresponds to the case without self-gravity {A4 = 10). 
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Fig. 5. Time evolution of (5E (solid line), PshSL (long-dashed 
line) and Sa (short-dashed line) for A^=10 and re = 1. 
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Fig. 6. Distributions of density perturbation Sp{z) at r = 0.2, 
0.4, 0.6 and 1.0. The abscissa is normalized by 
^0 = Cs/{\/2ttGpeM)- The ordinate is normalized by p^Ai^ ■ 



to be independent of the Mach number if we take f as di- 
mcnsionless time. Therefore, the growth rate of our calcu- 
lation is also expected to be independent of the Mach num- 
ber using the normalized time, t/tc, and the wavenunibcr, 

Next, we consider the transition epoch from an oscilla- 
tory to a growing mode. In figure 8, it is clearly seen that 
(5E changes from the oscillatory to the growing mode at 
the same value of rA^^/^. The filled circle indicates the 
transition epoch Tc^k predicted from the 1-zone model. In 
both of figures 8 (a)K = 1 and (b)K = 2, it is seen that 
ST, begins to grow at TgM^^^ ~ fc,K, where Tg corresponds 
to the inflection point of ST. We checked that Tc^k in 
the 1-zone model gives a correct transition epoch even in 
the model without the 1-zone approximation, at least for 
0.125^ k£ 8. 

5.3. Time-Evolving Dispersion Relations 

Dispersion relations at each instant of time are consid- 
ered. The growth rate at t is defined by 



(34) 



Figure 9a shows the dispersion relation derived from equa- 
tion (34) at each instant of time for Ai = 30. Figure 9b 
shows dispersion relations for a static layer correspond- 
ing to the density distribution at each instant of time (see 
GL65). When t < 1, the distinction between figure 9a 
and figure 9b is obvious. In the static layer, when r — > 0, 
(v— a^)niax and k:^}^^ are given by the free-fall growth rate 
and the thickness of the layer, respectively (see Appendix 
3). However, in the shocked layer, (V— i^^)max and fcmax 
converge to the zero as t 0. Thus, imposing shock 
boundary condition is important to understand the frag- 
mentation process in the shocked layer. Welter (1982) 
reported a similar property for a layer with shock on one 
side. 

5.4. Fragmentation Time: When Perturbation Becomes 
Non-linear 

Since the growth rate depends on both the time and the 
wavenumber (figure 9a; figure 2b), the above discussion is 
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Fig. 7. Cross sections of the perturbed layer for Ai = 10 and k = 2 at (a)r = 0.6 and (b)r = 1.0. The normalization Sz^h/^sh = 0.05 
at r = 0.6 is used. Contours denote the density perturbation normalized by p^Ai^, and arrows denote the velocity perturbation 
fields of which length represents their magnitude. The thick solid lines show Zgij = z^ij g + (52sh. The abscissa and the ordinate are 
normalized by zq as figure 6. The unperturbed shock position at r = 0.6 and 1.0 arc z^i^ g = 0.49^0 and 0.62zq, respectively. 




Fig. 8. The same figure as the figure 4 except that the abscissa denotes normalized time f = rAi which is defined in subsection 2.2. 
Filled circle shows the transition time -fc,fc which is predicted from 1-zonc model in subsection 2.2. 



not sufficient to predict the fragmentation scale. Thus, 
we consider the wavenumber that becomes non-linear at 
the earliest epoch. The time t^i when a mode becomes 
non-linear is defined by 

^S(Tni) _ 

TT-, — r = l- (35) 

The time Tni depends on the wavenumber and the initial 
amplitude of the perturbation. To focus on the evolu- 
tion after the transition epoch, Tc^k, we define the initial 
amplitude Sq at Tc^k., by 



Combining this definition with Eg oc t, we obtain 



-^0 



1. 



(36) 



(37) 



The arrival time t^i, which satisfies equation (37), is cal- 
culated as a result of our linear analysis. 



5.5. The Fastest Growing mode 

In this section, we determine the wavenumber that be- 
comes non-linear at the earliest epoch by using equa- 
tion (37). Figure 10 represents the dependence of Tni on 
the wavenumber for various Ai and initial amplitude, Sq. 
Figure 10 shows that the mode with k = 2 becomes non- 
linear at the earliest epoch for all M and Sq. Note that 
this wavenumber is different from one (k = 1, figure 2), 
which becomes unstable at the earliest epoch (see sub- 
section 2.2). For the fastest growing mode, n^i = 2, the 
dependence of Tni on the Mach numbers is shown in figure 
11. From figure 11, we can see that Tni (xAi^^^^ holds for 
each Sq. This Mach number dependence is simply under- 
stood by the fact that equation (33) docs not contain A4 
with normalized t/tc and fc/fc^. By fitting the dependence 
on 5q, t„i is finally given by 



T„i^2.4<So-°iX-i/2. 



(38) 



The error of this fitting formulae is less than 10% for 
10-3 < Jo < 0.1. 

Since the dependence on Sq is small in equation (38), if 
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(a) shocked layer 




(b) static layer 
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Fig. 9. (a) Dispersion relations derived from our calculation at each epoch for Ai = 30. (b) Dispersion relations of a static isothermal 
layer with the column density and an external pressure at each epoch (see GL65). The abscissa and the ordinate are normalized by 
the scale height Hq and 27rGpoOi respectively. 



(a)M=10 




So=0.1 — B- 

0.05 — ■- 
0.01 ■■■o- 



0.5 1 1.5 2 2.5 3 3.5 4 
K 



(b)M=30 




5g = 0.1 — B— 

0.05 — ■— 
0.01 ---o--- 



0.5 1 1.5 2 2.5 3 3.5 4 
K 



Fig. 10. Non-linear time defined by equation (37) for wavenumbcrs and initial amplitude for the case with A4 =(a)10, (b)30 



we neglect factors of order unity, Tni becomes ^ Ai ^1'^ . 
Therefore, i„i = T„i/V27rGpE ~ t|'°"'^/%/AT is inversely 
proportional to the growth rate at tc,K, where = 
l/V27rGpE- Note that tni is shorter than the free-fall 
timescale of the preshock region ~ t^""*^ and larger than 

the free-fall timescale of the layer ^^''^ ~ tf°'^'^IM. For 
larger Mach numbers, the difference among <ni, 
and becomes large. The origin of dependence oc 
1 /VTA is explained in subsection 2.2. In summary, k^i = 2 
and equation (38) are one of the main results in our linear 
analysis and that provide fragmentation length and time 
scales. 

6. Discussion 

6.1. Astrophysical Implications 

As an application of equation (38) , we seek criteria for 
collision-induced fragmentation. As a simplest case, we 
consider a spherical cloud with uniform density pE, the 
mass M and radius R. This cloud is specified by one 
parameter, ao, which is defined by 



— 1 1 1 — r 




—I I I i_ 



10 



100 



Mach number 



Fig. 11. Dependence of t,ii on Mach numbers for the case 
with K = 2. Initial amplitude (5o=0.01, 0.05, and 0.1 are de- 
fined in equation (37) 



5 c^R 



ao 



|Wo| 2GM' 



(39) 



Gravitational Instability of Shocked Interstellar Gas Layers 



9 



where Wq and IAq are the gravitational and thermal en- 
ergy of the cloud at the initial state, respectively. From 
equation (39) and M ~ 47rpE-R'^/3, the radius is given by 



2R: 



15c2 



27rG'pEao 



(40) 



We consider the collision of two physically identical clouds 
with velocity ±A^Cs relative to the center of mass. 

Fragmentation during cloud collision is expected if the 
following two conditions are satisfied (hereafter we adopt 
tni as fragmentation time). 

1. ^ Tconti The first condition is that fragmenta- 
tion time, Till, is shorter than the normalized colli- 
sion continuous time, Tconti, which is defined by 

Tconti = -rj- V^TtGPE ■ (41) 

From equation (40) and (41), wc obtain 



(42) 

Substituting equation (42) and our result of linear 
analysis [equation (38)] into Tni^Tconti, the first con- 
dition for fragmentation is written as 

ao<2Ml-^M~^. (43) 

2. 2Cs<conti ^ 2i? — A„i In order to assume the com- 
pressed region between two clouds as a layer, rar- 
efaction wave from the cloud surface must not erode 
in an unstable region with Ani, which is given by 



ISa, 



2ti 

Anl = -, — 211- 



, (44) 

fcni 2V27rG>EX 

This condition is expressed as 2csiconti Iz'^R^ \\\- 
With equation (40), (42) and (44), the above condi- 
tion is rewritten as 

2 



< 



(45) 



Therefore, two conditions, (43) and (45), provide fragmen- 
tation criteria on and M. 

Figure 12 represents the two conditions on the (ao,A^) 
plane. Fragmentation is expected in the dotted region. 
Figure 12 indicates that the case with ao — 1 can not be 
expected to fragment in all Mach numbers. This means 
that the collision of clouds in dynamical equilibrium does 
not induce fragmentation. By collision between clouds, 
only the case with small ag is expected to fragment the 
layer. In other words, only clouds with mass much larger 
than Jeans mass will fragment by cloud collision as long 
as isothermal. 

If cto ^ Ij each cloud collapses due to self- gravity. 
The collapse timescale is tdyn = ■\/37r/ (32Gpe), or Tdyn = 
•\/37r/4. The critical Mach number, A^cq, in which Tni 
equals to r^yn, is given by 

Meq~3.1,5o-°-2. (46) 



In Figure 12, thick dotted vertical line represents M = 
Mc(i- In the region with M ^ Mcq (ini ^ idyn), the dy- 
namical collapse of each cloud can not be ignored during 
fragmentation. Since our linear analysis docs not include 
collapse of each cloud, a non-linear multi-dimensional cal- 
culation is required for further an investigation. In the 
region with M ^ Mcq {tni ^ ^dyn), the layer is expected 
to fragment before an individual collapse. 

Nagasawa & Miyama (1987) investigated collisions be- 
tween self-gravitating isothermal equilibrium spheres sup- 
ported by an external pressure using three-dimensional 
calculations. They calculated cases by coalescence be- 
tween two clouds with ao ~ o(l) and ^ 10. As a result, 
they found that a single cloud forms and it becomes grav- 
itationally unstable. However, it does not fragment, but 
collapses as a single object. Our results are consistent 
with their results in ao ~ o(l) and M ^ 10. However, in 
order to compare our results and to discuss the fragmen- 
tation process of the layer which is formed by a collision 
between clouds with ao < 1, a much higher numerical res- 
olution than previous calculation is required. 

Boyd & Whitworth (2005) investigated the evolution 
of an oblate spheroidal cloud formed by a shocked layer 
using a semi-analytic method. Their clouds correspond to 
the unit wavelength in out model, while they included the 
effect of detailed thermal evolution. A cloud smaller than 
^ Cg/ y/2nGpE.M oscillate, and then collapse, while larger 
clouds collapse monotonically. Their results are consistent 
with ours. 

Finally, we comment on two possibilities of fragmenta- 
tion that is induced by cloud collision in an actual astro- 
physical environment. The first possibility is a collision 
between clouds with ao < 1 {M > Mj). From figure 12, 
the timescale of collision-induced fragmentation is smaller 
than that of the collapse of each cloud if collision velocity 
is ^ A^cq ~ 5 [equation (46)] . Therefore sheet fragmen- 
tation becomes dominant process in such a case of cloud 
collision. The second possibility is that ao becomes effec- 
tively smaller than order of unity after formation of the 
layer. This will be realized if the postshock region be- 
comes cooler than the preshock region by cooling, even in 
a collision between clouds with ao ^ 1. 

7. Summary 

We investigated the gravitational instability of isother- 
mal layers bounded by shock waves. The results of our 
investigation are summarized as follows. 

1. By imposing the shock boundary condition, the 
shocked layer is more stabilized than static layer, 
which is bounded by a constant external pressure. 
The stabilization influences the fragmentation pro- 
cess of the shocked layer. 

2. In our linear analysis, the epoch at which the per- 
turbations change from the oscillatory mode to the 
growing mode is well approximated by tc.K [see equa- 
tion (6)] derived from the 1-zone model. 

3. Our unsteady linear analysis can provide the correct 
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Fig. 12. Criterion of fragmentation. The solid and long- 
dashed lines correspond to r,ii = Tconti and 2cstconti = 2-R— Ani, 
respectively. The fragmentation region is shown to be 
in the dotted region. The short-dashed line represents 

Tnl = Tdyn ( = eq ) • The region with Tnl > Tdyn [M <Mcq) 

shows that the layer probably fragments during individual col- 
lapse. The region with t^i < r^yni-M > A^cq) shows that the 
layer can fragment before individual collapse. 

fastest growing mode, which becomes non-hnear at 
the carhcst epoch. The wavenumber of the fastest 
growing mode is kn\ = 2y/2TTGpEM/cs (khI = 2). 
This value, k^i = 2, is twice as large as the wavenum- 
ber. K = 1, which begins to grow at the earliest 
epoch. The time at which this mode becomes non- 
linear is given 

4. Collision-induced fragmentation is expected only 
when ao = 5c^R/2GM of the parent clouds is much 
smaller than 1 (figure 12). In the collision between 
clouds with Q!o < 1 and A4^5, the layer will frag- 
ment before individual collapse by self-gravity. 
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Appendix 1. The Shock Boundary Conditions 

A. 1.1. Unperturbed Shock Boundary Conditions 

In the rest frame of the postshock. wc derive jump 
conditions of isothermal shock without self-gravity. A 
schematic picture is represented in figure 13. The phys- 
ical variables in the preshock and the postshock regions 
are denoted by subscripts of "E" and "sh", respectively. 

To begin, jump conditions in the unperturbed state 
are considered. From mass and momentum conservation 
across the shock front and the isothermal condition, wc 
can derive the following equations: 



shock 
postshock region 
Psh 



front 
preshock region 

Pe 
-Ve 

< 

dzsh /dt 



Zsh 



z axis 



Fig. 13. Schematic picture of shock front in the rest frame 
of postshock. 

Psh + Psh{-Zshf = Pe + Pe{~Ve - ish)', (A2) 
and 

Psh Pe 

We define the flow velocity V^h in the rest frame of the 
shock front by Vsh = Ve + ish- Hence, by solving equations 
(Al), (A2) and (A3), wc obtain 

Psh = PeM^, 

Zsh 



(A3) 



(A4) 



M' 
and 

Ve^IM 



M 



A. 1.2. Perturbed Shock Boundary Conditions 



(A5) 
(A6) 



Linearized jump conditions are derived. The posi- 
tion of the shock front is defined by Zsh(x,t) = Zsh.o(^) + 
5zsh(i)e*'"^. Linearized mass and momentum conservation 
across the shock front are given by 

Dpsh{-Zsh) + PshfliDv^sh - Szsh) = -pESzsh, (A7) 
and 

DPsh{cl + igh) - 2/9sh,0-2sh(-Dw^,sh - SZsh) 

= 2peV,^5z,^, (A8) 

where operator "Z?" represents Lagrangian displacement. 
From equations (A7) and (A8), we have 

M d5zsh 



Dpsh = 2pE- 
and 

£'Wz.sh = ( 1 



dt ' 

1 N dSzsh 
M^j dt ■ 



(A9) 



(AID) 



Psh(-ish) = Pe{-Ve - ish), 



(Al) 



The relationship between the Lagrangian and Eulcrian 
displacement is given by 
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Dpsh^Dpsi,{zsh,o + 6zsh)^6psh + Szsh — — 



dpo 
dz 



.(All) 



Next, we consider the velocity perturbation in the x- 
direction 5vx,sh- Momentum conservation for the tangen- 
tial direction of the shock front gives 

d5zsh 



dx 



(A12) 



We adopt the same boundary condition of the gravita- 
tional potential perturbation as in Welter (1982), 



dz 



kScjjsh + 47rG(psh - Pe) = 0. 



(A13) 



Appendix 2. Dispersion Relation of a Shocked 
Layer under 1-zone Model 

The averaged physical quantities arc defined by 



pdz and = 



pvxdz. 



(A14) 

The Oth-ordcr equation is given by 

atEo-pEl/sh, (A15) 

where the subscript "0" represents the unperturbed vari- 
able and differential operators d/dt and d/dx are abbre- 
viated by dt and dx, respectively. Let us consider the 
perturbation variables 6Jx and 5L, which are defined 

by 



(A16) 
(A17) 

2z,i,{x,t). (A18) 

The column density perturbation, JS, can be divided into 
two contributions, as follows: 



S(.T,<) = I]o(t) + (5I](a;,t), 
Jx{x,t) = SJx{x,t), 
and 

L{x,t) = La{t) + dL{x,t), L{x,t) 



5Ti = PshSL + 6a, where Sa 



Spdz. 



(A19) 



^sh,0 



The linearlized equation of mass conservation is given 

by 

dtSE = pEdtSL - dxSJx ~ -dxSJx. (A20) 

Since ST, ^ Psh^L ^ pe^L, linearlized equation of momen- 
tum conservation in the x-direction is given by 

dt6Jx = PEVshi~VEdx6L) - c%Sa + 27rGTodxSJ:/k 
- -PeM^OxSL - cldxSa + 2T:GTodx5T/k. 
= -cldx5T. + 2TiGT.Qdx5T./k. (A21) 

The second term in the first line of the equation (A21) rep- 
resents the momentum flux as a result of the shock bound- 
ary condition. The second line of the equation (A21) is 
derived using Aicg. The third line is derived by using 
Psh = Pe-M.^ and equation (A19). 

Neglecting the time-dependence of So(^)i from equation 
(A20) and (A21), the approximate dispersion relation is 
derived as 



27rGfcI]o(i)- 



(A22) 



Appendix 3. A Layer 
External Pressure 



Confined by 
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Constant 



We investigate a highly compressed static layer that is 
mainly confined by constant external pressure, Pb- In this 
situation, the layer has almost uniform density, denoted by 
Po- The momentum conservation equation in x-direction 
is given by 



dtSJx 



%Sa + 2TrGTodxST/k, 



(A23) 



where the second term of first line of equation (A21) is 
absent due to the constant pressure boundary conditions. 
Because the layer is very thin, Sp is approximately deter- 
mined only by boundary value of Sp as Sa ~ Sp{z = Zh)LQ, 
where Zb is the position at the boundary of the layer. The 
constant pressure boundary condition is given by 

Dp = <5p(zb,o) + Szk{d,po),^„ = 0, (A24) 

(see GL65). Because the unperturbed state is hydrostatic 
equilibrium, (9z/3o)zb = -Po{dz(l)o)zi,/c^ = -2ttGpILo/c^. 
Therefore, the ratio of Sa to poSL is written as 

Sa/poSL - {Lo/Hof < 1. (A25) 

From SY, ^ poSL and equation (A25), the second term of 
the equation (A23) becomes 

2 



-ddxSa^-ctij^^] dxSj: 



L 



(A26) 



5o = 2^Gpooio-27rGSo, (A27) 



Thus, the pressure term is ^ (Lq/Hq)'^ times as small 
as that of the shocked layer [equation (3)]. Moreover, 
Cs{L()/Hq) is rewritten with equation (7) as 

h. 

Ho 

where Sq — Poo^o- This corresponds to the phase velocity 
of gravity wave in shallow water. Finally, the dispersion 
relation is given by 

2 

27rGfcSo - 27rGEo(Lofc^ - fc)(A28) 



OJ 



Ho 



2, 2 



From equation (A28), (v— Ci^)max is almost the free-fall 
growth rate of the layer and Amax is nearly equal to the 
thickness ~ Lq. The smaller scale perturbation can grow 
much faster than in the shocked layer because the pressure 
effect is small. 

From the above simple analysis, we understand the fol- 
lowing things. For (Lo/iJo)~l, the difference of boundary 
condition does not influence on the evolution of the per- 
turbations, and the dispersion relations become the same 
form as equation (A22). However, for (Lq/Hq)^ <^ 1, the 
constant-pressure boundary condition is in completely dif- 
ferent from that of the shock boundary condition. 

Appendix 4. Validity of Approximation (t>z,o — 0) 

The validity of the approximation Vzfl = is considered 
in this section. The perturbation equations containing 
Vzfi are given by 
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Fig. 14. (a): Analytic expression (A32) is shown for A4 ■ 
(A32){vzfi), equation (3)(tj) and equation 



-A, 10, 30, 50. (b): Plots of v^fi / (zgi^uj). Each variable is from equation 



at az 



(A29) 



d d 



^ -ikpov^fiSvj;- Sp-^ - po-g^, (A30) 



and 



d d 

-g^ipoSvx) + -g^ipoVz.oSvx) = ~ikcl5p~ikpo6(t>. (A31) 

For \z\ < |zsh|, the velocity l^z^ol-z)! is smaller than 
bx.oC^sh)! at shock front. If we assume that the shock 
front position is given by equation (11), we can express 
Vz,oizsh) approximately from the jump conditions 



M \l-fiT)/M^ 
where / is given by 



+ ./W-1 



(A32) 



/(r) :=:TCOsh-Vl + rV(l + T)3/2 



(A33) 

Figure 14a shows the time evolution of |uz,o(^sh)| for var- 
ious Mach numbers. This figure indicates that the maxi- 
mum value of |wz.o(-2sh)| is slightly less than, or comparable 
to, ^ 0.3cs/A^. The average velocity, |uz,o|, of the layer is 
clearly smaller than this value. 

Let us consider the contribution of Vzflizsh) for the 
perturbation equations. To begin with, we compare the 
two terms Spvz.o and poSvz in perturbed momentum in 
the z-direction. The ratio of the two terms is given by 
iSpvz.o)/ipSvz,o) < 0.3cs{Sp/Svzfi)/M^ < 1; thus, we can 
approximate the perturbed momentum in the z-direction 
as pqSvz- Therefore, we can neglect the third term of 
equation (A29) and the second term of equation (A30). 
Moreover, we can neglect the fifth term of equation (A30), 
which is smaller than the third term by 0.3/ A^^. In equa- 
tion (A30), the ratio of the fourth to the first terms is 
given by Vz,o/{Lqlo). This ratio is small value if Mach 
number is large (see figure 14b). For the same reason, we 



can neglect the second term of equation (A31) compared 
with the first term. The ratio of the sixth term of equation 
(A30) to the first term is given by {k/uj)vz,o{Svx/5vz) ^ 
{0.3/M){6vx/Svz)- This value is also small if the Mach 
number is large. From the above discussion, Vz,o can be 
assume to be zero. 
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